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ABSTRACT 

The universal baryonic mass fraction can be sensitively constrained using 

X-ray observations of galaxy clusters. In this paper, we compare the baryonic mass 
fraction inferred from measurements of the cosmic microwave background with the gas 
mass fractions (/gas) of a large sample of clusters taken from the recent literature. In 
systems cooler than 4 keV, /gas declines as the system temperature decreases. However, 
in higher temperature systems, /gasC^^soo) converges to « (0.12 ± 0.02) (ft,/0.72)~^'^, 
where the uncertainty reflects the systematic variations between clusters at rsoo • This 
is significantly lower than the maximum-likelihood value of the baryon fraction from 
the recently released WMAP 3-year results. We investigate possible reasons for this 
discrepancy, including the effects of radiative cooling and non-gravitational heating, 
and conclude that the most likely solution is that is higher than the best-fit WMAP 
value (we find fi™ = O.SOlg Jg), but consistent at the 2a level. Degeneracies within 
the WMAP data require that cts must also be greater than the maximum likelihood 
value for consistency between the data sets. 
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1 INTRODUCTION 

For over a decade now cluster gas mass fractions as inferred 
from X-ray observations have been used as a probe of the 
universal ratio of baryon to total matter densities, S7i,/f2,„ 
(e.g.. White et al. 1993; David et al. 1995; Evrard 1997; 
Mohr et al. 1999; Roussel et al. 2000; Allen et al. 2002; Lin 
et al. 2003; Ettori 2003; Allen et al. 2004). Supplementing 
these gas mass fractions with constraints on SI;, from, e.g., 
cosmic microwave background (CMB) measurements or a 
combination of Big Bang Nucleosynthesis (BBN) predictions 
and D/H measurements from high redshift quasars, there- 
fore allows one to measure the total matter density Q.m- The 
reliability of this test rests on the assumption that clusters 
have been able retain the original baryon inventory assigned 
to them in the early universe. So-called "non-radiative" cos- 
mological simulations, which include a hydrodynamic treat- 
ment of the baryons but neglect sources or sinks such as 
radiative cooling, star formation, and feedback, indeed in- 
dicate that clusters retain nearly all their baryons until the 
present day (e.g., Frenk et al. 1999; Kay et al. 2004; Grain 
et al. 2006) . The same is generally true for simulations with 



E-mail: i.g.mccarthy@durham.ac.uk (IGM) 



cooling and feedback. Although the fraction of baryons in 
the hot phase depends strongly on the model, most recent 
simulations predict a mild increase in the hot gas fraction 
with cluster mass, and little evolution with redshift (e.g., 
Tornatore et al. 2003; Kravtsov et al. 2005; Ettori et al. 
2006). 

Although the cluster baryon fraction test has been ex- 
amined previously in many studies, there are several good 
reasons for revisiting it. First, new high-quality data ob- 
tained from the Chandra and XMM-Newton telescopes now 
allow us to probe both the surface brightness and tempera- 
ture profiles of clusters out to large radii. As a result, both 
the statistical and systematic observational uncertainties on 
the gas mass fraction are substantially improved. Second, 
much improved (e.g., K-band) measurements of the stellar 
content of clusters are now available. Third, cosmological 
simulations can now robustly predict the baryon fractions 
within rsoo, which is roughly the same radius the latest X- 
ray measurements reliably extend out to. Fourth, analyses 
of mock observations of realistic cosmological simulations al- 
low us to, e.g., quantify the observational bias introduced by 
assuming strict hydrostatic equilibrium (HSE) in the deriva- 
tion of X-ray gas mass fractions. Finally, analysis of the re- 
cently released WMAP 3-year cosmic microwave background 
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(CMB) data has yielded very tight constraints on the univer- 
sal baryon fraction (Spergel et al. 2007). As a result, there 
is now an excellent opportunity to check for consistency (or 
lack of) between these two very different tests. We note that 
Ettori (2003) previously found evidence for a slight discrep- 
ancy between the universal baryon fraction inferred from the 
first year WMAP data and cluster baryon fractions inferred 
from BeppoSAX X-ray data. It is interesting to see whether 
or not this discrepancy remains when one makes use of the 
better observational data and improved theoretical predic- 
tions that are now available. 

In this paper we analyse the best literature X-ray data 
in a homogeneous way, and compare the results with the 
universal baryon fraction inferred from the recently released 
WMAP 3-year data (Spergel et al. 2007). 

Unless otherwise stated, we assume a ACDM cosmology 
with h = 0.72. 



2 CLUSTER GAS FRACTIONS 

We select high quality data from a few recent studies in the 
literature where the mass profiles have been reliably esti- 
mated out to ~ rsoo. Vikhlinin et al. (2006), hereafter V06, 
have measured the gas and total mass profiles for a sample 13 
relaxed "cool core" observed with Chandra. From this sam- 
ple, we select all but two cluster^. In addition to the above, 
we also select a sample of 10 relaxed "cool core" clusters 
observed with XMM-Newton for which Pointecouteau et al. 
(2005), Arnaud et al. (2005) and Pratt et al. (2006a) (here- 
after collectively referred to as PAP) have computed gas and 
total mass profiles. Therefore, in total we have compiled a 
sample of 21 sets of mass profiles from 19 different clusters 
(i.e., the samples have two clusters in common). 

Both V06 and PAP computed the gas and total mass 
profiles of their clusters in a similar manner, by fitting para- 
metric forms of the gas density and (3D) temperature pro- 
files of the ICM (see eqns. 3 and 6 of V06 and Appendix 
A of Pratt et al. 2006a) to the observed, projected surface 
brightness and temperature profiles. The purpose of these 
parametric models is to produce a smooth description of 
the data and hence reduce the noise in the spatial deriva- 
tives; the exact form of the models is irrelevant as long as 
they fit the data well. Note that the assumption of smooth 
gas in HSE can lead to a small overestimate of the baryon 
fraction (Mathiesen et al. 1999; Mohr et al. 1999; Rasia et 
al. 2006), which strengthens our conclusions. 

We use the parametric models and associated best-fit 
parameters listed in V06 and PAP to reconstruct the ob- 
served mass profiles. In both studies the total mass distri- 
butions were fitted with the Navarro, Frenk & White (1997, 
NFW) profile derived from cosmological dark matter simu- 
lations. Both V06 and PAP demonstrate that the NFW pro- 
file fits their data exceptionally well, with an inferred mass- 
concentration parameter relation that is in good agreement 
with the results of cosmological simulations. 

^ We exclude A2390, which V06 have demonstrated to be highly 
asymmetric owing to the presence of a set of large X-ray cavi- 
ties nearly 400 kpc in diameter, and the low-temperature system 
USGC S152, for which V06 do not provide enough information to 
reconstruct its mass profile. 
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Figure 1. Comparison of the observed integrated gas mass frac- 
tions [i.e., /gas = Mgas(< A)/Mtot(< A)] as a function of over- 
density (A) with WMAP 3-year constraints on the universal 
baryon fraction Qi,/Qm (assuming a flat power-law ACDM cos- 
mology). Solid and dashed lines represent fits to the Chandra data 
of V06 and XMM-Newton data of PAP, respectively. Red and 
black lines represent systems with mean spectral temperatures 
below and above 4.5 keV, respectively. The magenta line repre- 
sents PAP's fit to A1413. Also shown (data points), are weighted 
averages of the gas mass fraction determinations of Ettori et al. 
(2002) (blue squares), Allen et al. (2004) (green triangle), Sadat 
et al. (2005) (magenta circle), LaRoque et al. (2006) (yellow pen- 
tagon), and Afshordi et al. (2007) (open cyan circle). Typically, 
the associated statistical measurement uncertainty for individual 
clusters is of order 10% for overdensities within which the tem- 
perature and surface brightness profiles can be reliably measured. 
The vertical dotted line indicates an overdensity of 500 (i.e., cor- 
responds to rsoo). The shaded cyan region corresponds to the 
68% confidence region for Qf,/Qm from WMAP, while the long 
dashed line shows the best fit value (Spergel et al. 2007). 



Presented in Figure 1 is a comparison of the observed 
integrated gas mass fractions as a function of overdensity, 
A, where A = 3Aftot(< r)/ [47rr^pcrit(2)] . In general, the 
clusters all show a mildly rising gas fraction with decreas- 
ing overdensity. However, there is considerable scatter in 
the gas fraction at fixed overdensity that is worth exploring. 
First, it is evident that the gas mass fractions measured by 
PAP are systematically larger than those measured by V06. 
Given that the gas density profiles from the two samples are 
quite similar (see McCarthy et al. 2007), the implication is 
that the temperature profiles measured by PAP and V06 
differ systematically from each other. Indeed, V06 generally 
find temperature profiles that decline relatively rapidly with 
radius, dropping by roughly a factor of 2 from the peak 
(at r ~ 0.1 — 0.2r2oo) to 0.5r2oo (see also Vikhlinin et al. 
2005), while PAP find a much more gradual decline, with 
some clusters showing approximate isothermality. Through 
the equation of HSE, a flatter temperature gradient trans- 
lates into a reduced normalisation of the total mass profile 
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Figure 2. Comparison of the observed integrated gas mass frac- 
tions within an overdensity of A = 1500 as a function of mean 
ICM temperature with WMAP 3-year constraints on the univer- 
sal baryon fraction Qi,/Qrn- Solid squares represent the Chandra 
data of V06 and open triangles represent the XMM-Newton data 
of PAP. 

and hence an increased gas mass fraction. While it would be 
useful to sort out the exact nature of the temperature dis- 
crepancj0, we point out that the typical level of difference 
between the two is relatively small (compared to the offset 
of both from the WMAP result, see §3). The one exception 
to this appears to be A1413, which was observed by both 
V06 and PAP. Within rsoo, V06 find a total mass that is 
roughly 50% larger than that found by PAP. Increasing the 
total mass of A1413 by this factor would bring it more in 
line with the other systems studied by PAP. However, we 
note that good agreement between PAP and V06 is found 
for A1991, the only other system in common between the 
two samples. 

Fig. 1 also shows that the gas fraction depends on sys- 
tem temperature, which presumably reflects total system 
mass. In particular, systems with mean temperatures be- 
low 4.5 keV (red lines) have systematically lower values for 
/gas within virtually all overdensities compared with hot- 
ter systems (black lines). This trend has been noticed pre- 
viously (e.g., Mohr et al. 1999; Roussel et al. 2000; Neu- 
mann & Arnaud 2001) and would seem to indicate that 
non-gravitational heat sources are relatively more impor- 
tant in mediating the properties of the ICM in low mass 
than in high mass clusters. In Figure 2, we plot the gas 

^ It now appears that the discrepancy between temperature pro- 
files derived with Chandra and those with XMM-Newton may be 
nearly resolved. Using a larger sample of clusters observed with 
XMM-Newton and an improved model for background subtrac- 
tion, Pratt et al. (2006b) now find steeper temperature declines at 
large radii which are quite similar to those of V06 (G. W. Pratt, 
private communication). See also Vikhlinin et al. (2005). 



mass fraction within an overdensity of 1500 (typically cor- 
responding to a physical radius of ~ 800 kpc) versus mean 
system temperature. We have chosen this radius so that all 
the clusters in the sample (including the coolest systems) 
have reliable data. Note that if non-gravitational heating 
is unimportant at risoo it will be even less significant at 
rsoo. A trend with temperature is clearly visible at low tem- 
peratures. However, above ~ 4 keV both the Chandra and 
XMM-Newton data show no evidence for a further increase, 
and the gas fraction appears to have converged, suggest- 
ing that non-gravitational physics is largely unimportant for 
the most massive systems (at least out at large radii). This 
is fully consistent with previous results based on ROSAT 
and ASCA data (Roussel et al. 2000; Lin et al. 2003), as 
well as newer XMM-Newton data (Sadat et al. 2005). It is 
also consistent with the findings of McCarthy et al. (2007), 
that above ~ 3 keV the gas density and entropy profiles 
of "cool core" systems follow the scaling expected for self- 
similar gravitational heating at large radius. 

We have elected to focus on the data of V06 and PAP 
since their measurements extend out to r-5oo. As noted 
earlier, cosmological simulations can now reliably predict 
the baryon fraction within this radius. However, measur- 
ing the properties of the ICM at such large radii is not easy. 
The rapidly declining surface brightnesses of the clusters at 
large radii means that an extremely careful treatment of the 
background is crucial. For this reason, a number of groups 
have limited their analyses to smaller radii (higher surface 
brightnesses) in an attempt to mitigate any potential sys- 
tematic effects associated with the background modelling. 
While we are encouraged by the similarity of the results of 
V06 and PAP at large radii (note also that Chandra and 
XMM-Newton have different background characteristics), it 
is useful to check these results against other (mainly small- 
radii) data from the recent literature. 

With the above in mind, we include some other re- 
cent X-ray gas mass fraction determinations (within vari- 
ous fixed overdensities) of hot clusters from the literature 
(data points) in Figure 1. In particular: the (yellow) pen- 
tagon represents a weighted average of /gas (''2500) for the 38 
clusters (spanning 0.14 < z < 0.89) observed by LaRoque et 
al. (2006) with Chandra [note that LaRoque et al. actually 
use both Chandra X-ray data and BIMA/OVRO Sunyaev- 
Zeldovich effect data to constrain /gas(r25oo)]; the (green) 
triangle represents a weighted average of /gas(r250o) for the 
26 clusters (spanning 0.07 < z < 0.9) observed by Allen 
et al. (2004) with Chandra; the (blue) squares represent 
weighted averages of /gas(r250o) and /gas (r 1000) for the 22 
clusters (spanning 0.01 < z < 0.1) observed by Ettori et 
al. (2002) with BeppoSAX, and; the (magenta) circle repre- 
sents a weighted average of /gas (rsoo) for the 8 clusters (at 
z^O.5) observed by Sadat et al. (2005) with XMM-Newton. 

Although there is some scatter in the averaged /gas val- 
ues from these different studies. Fig. 1 demonstrates that 
they are, by and large, consistent with the profiles of V06 
and PAP. The one exception is the data of Ettori et al. 
(2002), which are systematically lower than the others. Since 
this is the data Ettori (2003) used to show there was only 
a very small discrepancy between cluster baryon fractions 
and the first-year WMAP results, we discuss the studies of 
Ettori et al. (2002) and Ettori (2003) in a bit more detail in 
§3. 
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X-ray observations are not the only way one can con- 
strain the gas mass fractions of clusters. Increasingly, the 
Sunyaev-Zeldovich (SZ) effect is also being used for this pur- 
pose. For example, LaRoque et al. (2006) found very similar 
results for /gas(?'250o) when fitting to both Chandra X-ray 
data and BIMA/OVRO SZ effect data simultaneously (yel- 
low pentagon in Fig. 1) or when fitting to the SZ effect data 
alone. More recently, Afshordi et al. (2007) have stacked the 
SZ effect signals of 193 clusters in the WMAP 3-year data to 
constrain the thermal energy of the ICM. With the aid of hy- 
drodynamic simulations, they converted this thermal energy 
into an estimate of the gas fraction of clusters within 7-200 • 
Their result, represented by the open (cyan) square in Fig. 
1, is in excellent accordance with the X-ray-derived results. 
This is very encouraging since the systematics involved in 
the derivation of this result are obviously very different from 
those associated with X-ray measurements. 



3 IMPLICATIONS 

In cosmological simulations the baryon fraction within rsoo 
has converged to 90 — 95% of the universal value in the case 
of non-radiative simulations (e.g., Frenk et al. 1999; Kay et 
al. 2004; Grain et al. 2006) and slightly higher than this 
when radiative cooling and star formation is included (e.g., 
Kravtsov et al. 2005; Ettori et al. 2006). It is only possible 
to drive a substantial fraction of the baryons beyond rsoo 
if the energy input from non-gravitational heating (such as 
AGN powered jets and bubbles; e.g., Ghurazov et al. 2002) 
is comparable to the binding energy of the cluster gas (we 
quantify the energy requirement in §3.2). In the absence of 
such high levels non-gravitational heating, the baryon frac- 
tions derived at rsoo should be nearly representative of the 
universal value. 

Independent measurements of the universal baryon frac- 
tion can be derived from the power spectrum of 
CMB anisotropics. In Figs. 1 and 2 we show the value of 
0.bl^m constrained by the WMAP 3-year data (Spergel et 
al. 2007). From Fig. 1 it is immediately apparent that, with 
the exception of PAP's estimate of /gas for A1413, the obser- 
vational gas mass estimates do not achieve consistency with 
the WMAP measurements of the baryon fraction within any 
(observable) overdensity. This is true despite the fact that 
the XMM-Newton data of PAP extend out to nearly rsoo 
(note that, typically, r500 '-^ 1 Mpc) and slightly beyond 
for the Chandra data of V06. At the lowest observable over- 
density (largest radii), the data indicate a gas mass fraction 
that is roughly 40% lower than Qi/^lrn inferred from WMAP 
3-year data. Fig. 2 demonstrates that the observed gas mass 
fractions at > risoo have converged for systems above ~ 4 
keV, and still lie well below the WMAP 3-year constraints. 
We now examine the possible origins of this sizeable dis- 
crepancy. For ease of discussion, we break up the possible 
solutions into three broad categories. 

3.1 Stars and cool baryons 

X-ray data by itself constrains only the fraction of a cluster's 
mass in the form of hot [T > 10®'^ K) gas. A proper com- 
parison to the universal WMAP baryon fraction therefore 
requires that we take into account the fraction of cluster's 



baryons locked up in cool gas (that doesn't emit X-rays), 
stars and baryonic dark matter. If clusters manage to sig- 
nificantly cool « 40% of their baryons this would poten- 
tially resolve the discrepancy described above. Indeed, cos- 
mological simulations that take into account the effects of 
radiative cooling and star formation demonstrate that clus- 
ters can potentially cool out such large quantities of their 
baryons (e.g., Dave et al. 2002; Kravtsov et al. 2005). But 
such simulations are at odds with near-infrared (K-band) 
observations of clusters, which typically indicate that the 
total (resolved) stellar mass is at most 5-14% of the IGM 
gas mass in hot (> 4 keV) clusters and only slightly higher 
than this in cooler system^ (e-g-, Balogh et al. 2001; Lin et 
al. 2003). Recent deep optical observations limit the contri- 
bution of diffuse intracluster light to between ~ 10 — 30% 
of the total stellar luminosity (e.g., Gonzalez et al. 2005; 
Zibetti et al. 2005). Low mass stars and brown dwarfs (or 
"rocks") are also likely to be a significant, undetected com- 
ponent of the mass budget, but there is no strong evidence 
that they are more abundant than expected from standard 
initial mass functions used to extrapolate the observed stel- 
lar mass functions (Fuchs, Jahreiss & Flynn 1998; Gizis et al. 
2000; Lucas et al. 2005; Levine et al. 2006). Radio, infrared, 
and ultraviolet surveys for atomic, molecular, and ionised 
gas (respectively) limit the contribution of cool (T < 10^'^ 
K) diffuse baryons to less than a percent or so of the hot 
X-ray-emitting IGM (e.g., O'Dea et al. 1998; Donahue et al. 
2000; Edge 2001; Edge & Frayer 2003; Bregman et al. 2006). 
Finally, Ettori (2003) have suggested that a significant frac- 
tion of cluster baryons could lie at temperatures between 
10^ < T < 10^ K and they cite possible evidence for this in 
the form the observed "soft X-ray excess" . However, it now 
appears that this excess was the result of inadequate mod- 
elling of the Galactic foreground (Bregman & Lloyd-Davies 
2006). Furthermore, it is unclear physically how such large 
quantities of gas in this temperature range could resist cool- 
ing down to 10'' K or being mixed or heated to the ambient 
IGM temperature. 

Furthermore, one can place limits on any hidden cold 
baryonic component by considering its gravitational effects 
on the properties of the IGM. In particular, if the hidden 
material is centrally concentrated, as one would expect if 
its origin is linked to cooling, this will deepen the poten- 
tial well which, in turn, will heat the IGM through simple 
gravitational compression. The observed temperature pro- 
files can therefore be used to place constraints on the amount 
of hidden cold material in clusters. We have tested this as 
follows. In particular, we assume the dark matter follows a 
NEW distribution with a typical concentration parameter 
C200 = 4. We further assume the gas density profile follows 
a distribution typical of observed "cool core" clusters (see 
V06; McCarthy et al. 2007). Finally, we try different hid- 
den cold mass distributions. The predicted IGM tempera- 
ture distributions are computed by placing the gas in HSE. 
In the case where there is no hidden component, we find 

^ We note that some previous studies (e.g., Ettori 2003; Voevod- 
kin & Vikhlinin 2004) estimated slightly higher stellar fractions 
(typically 15-20%) based on V- and B-band cluster observations. 
However, due to uncertain age and metallicity effects, conversion 
to stellar mass from the V- and B-bands is less reliable than con- 
verting from the K-band. 
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the resulting temperature profiles quantitatively match the 
observed ones. This is not unexpected, given the results of 
McCarthy et al. (2007) (see §3.1 of that study). However, 
the predicted temperature profiles are quite sensitive to the 
addition of a centrally-concentrated cold component. In par- 
ticular, if we require that the predicted temperature profiles 
are within the observational uncertainties (typically, binned 
temperature profiles have a 10% measurement uncertainty) , 
we can rule out a cold component weighing > « 3 x 10^^ Mq 
within the central 50 kpc or > ~ 5 x 10^^ Mq within the 
central 100 kpc. (Note that our results at a fixed radius are 
not sensitive to the distribution of the cold component, only 
to the total enclosed mass within that radius.) This is small 
compared to the ~ 5 x 10^^ Af© required to resolve the clus- 
ter vs. CMB discrepancy for a massive lO^'^ Mq cluster. The 
only way to reconcile the cluster vs. CMB discrepancy and 
maintain consistency with the observed temperature pro- 
files is if the hidden component is very spatially-extended. 
However, such a configuration seems physically contrived. 

It therefore appears that the stellar/cool baryon contri- 
bution to the total cluster baryon budget too small to ac- 
count for the cluster vs. CMB discrepancy. In particular, the 
resolved stellar/cold baryon component is a approximately 
factor of 3 smaller than required. Moreover, we have demon- 
strated that any hidden component cannot be very massive, 
otherwise this would violate the observed temperature pro- 
files of relaxed clusters. 



3.2 Non-gravitational heating 

As discussed above, the small fraction of cooled baryons 
observed, relative to the predictions of cooling-only simu- 
lations, implies that some form of non-gravitational heat- 
ing ("feedback") is at work in the ICM. Without significant 
feedback to prevent this overcooling, theoretical models are 
unable to account for the observed cut-off of the galaxy lumi- 
nosity function at the bright end or the fact that BCGs are, 
by and large, "red and dead" (e.g., Benson et al. 2003; Bower 
et al. 2006), nor would it be possible to explain the lack of X- 
ray emission from intracluster gas with temperatures below 
about 1 keV (e.g., Peterson et al. 2003). Non-gravitational 
heating also appears to be necessary to account for the X- 
ray scaling properties of clusters (e.g.. Kaiser 1991; Evrard & 
Henry 1991; Babul et al. 2002; Voit et al. 2002; McCarthy et 
al. 2004) . Injecting thermal energy into the ICM will cause it 
to expand (while leaving the dark matter in-tact) and there- 
fore will reduce the gas mass fraction within a given radius. 
Can this heating explain the discrepancy between cluster 
and CMB measurements? 

To test this, we have computed the bulk energy required 
to transform clusters with the universal baryon fraction (all 
in hot gas) into the observed systems. For our baseline model 
clusters, we assume the gas traces the total matter at all 
radii and, therefore, within any radius the integrated gas 
mass fraction is always the best-fit WMAP 3-year value 
fit/f^m = 0.176. The total mass profiles (which are dom- 
inated by dark matter) are assumed to be the same as those 
measured by V06 and PAP, thus we construct a baseline 
model cluster for each of the observed systems. The tem- 
perature profiles are determined by placing the gas in HSE. 
Calculating the total energy of the gas in these systems (i.e., 
the summation of the total internal and potential energies) 



is then straightforward. For the observed systems, we extrap- 
olate the gas and total density profiles beyond the maximum 
radius to which they can be observed, until the integrated 
gas mass fraction is the universal ratio. (Note that we also 
extrapolate the baseline models to this radius to ensure that 
both the baseline and observed systems have the same inte- 
grated gas and total masses.). We assume the total density 
profiles continue to follow the NFW form fit by V06 and 
PAP. We try various different powerlaw extrapolations for 
the gas density profiles with plausible indices ranging from 
to —2.5. We find that the minimum energy required to 
convert the baseline models into the observed systems oc- 
curs when the gas density is constant with radius outside 
the maximum observable radius. This configuration is per- 
haps unlikely, but it does provide a useful lower limit to the 
amount of heating required. Like the baseline models, we 
place the gas in the observed systems in HSE. Reassuringly, 
we verify that the resulting temperature profiles are in good 
agreement with the observed profiles. The total (minimum) 
energy required to heat the ICM is just the difference of the 
total energy of the observed and baseline model systems. 

For hot (Tspec > 4 keV) clusters, we find that a sub- 
stantial amount of energy is required, ranging between 
6 — 45 X lO''^ ergs with a mean of ^ 2.2 x lO''^ ergs (as- 
suming pgaa is constant outside the maximum observable 
radius - i.e., the minimum required energy). It is interesting 
to compare this minimum energy estimate with the energy 
that can potentially be deposited by AGN, the most power- 
ful source of non-gravitational heating we know of in clusters 
of galaxies. (We note that the energy available from AGN 
far exceeds that available from stellar feedback.) We esti- 
mate the amount of AGN energy available to be tapped as 
follows. First, we convert Lin et al. (2003)'s observed rela- 
tionship between stellar mass fraction and total mass within 
''500 (see equation 10 of that studjQ) into a stellar — total 
mass relation [i.e., Mstar(T'5oo) — Afsoo] assuming a gas mass 
fraction of 0.12 within rsoo (see Fig. 1). This relation can 
be converted into a relationship between the total mass in 
black holes within rsoo by (optimistically) assuming that the 
entire stellar mass is contained in bulges and adopting the 
black hole — bulge mass (A/bh — A^buigc) relation of Haring 
& Rix (2004). This leads to the following AfBH(r5oo) - Afsoo 
relation: 



logio [A/bh (rsoo) (^©)1 = 1-12 logi 



+ 0.84 logi 



ATs. 



5 X 1014 A/q 



/gas (rsoo) 
0.12 



+ 10.229 



(1) 



Finally, this is converted into an estimate of amount of 
AGN energy available via -Eagn = eAifBH(rsoo)c^. 

In Figure 3, we present a comparison of the minimum 
specific energy required to resolve the cluster vs. CMB dis- 
crepancy with the energy available to be tapped in black 
holes. (To facilitate comparison with the observations, we 
calculate the specific energy by dividing the total required 



* We have slightly adjusted this relation by scaling up the total 
masses of Lin et al. by 1.26 to account for the normalisation differ- 
ence between the ASCA total mass-temperature relation assumed 
by Lin et al. and the more accurate Chandra relation measured 
by V06. 
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of magnitude short of the required minimum to reduce a 
massive cluster's baryon fraction from the universal WMAP 
value to the observed fraction. Therefore, even if a typical 
cluster experiences 10 such powerful outbursts over its life- 
time (say, once per Gyr over 10 Gyr) the energy injected into 
the ICM still falls short of the minimum required energy by 
up to an order of magnitude. 

We therefore conclude that AGN heating is a highly im- 
plausible, but perhaps not impossible, solution to the cluster 
vs. CMB discrepancy. In addition to the exceptionally large 
energy requirements, we point out that the heating must be 
distributed in just such a way as to explain the convergence 
trend in Fig. 2 and the fact that the ICM proj)(>rtics at large 
radii in massive clusters follow the gravitational self-similar 
scalings (McCarthy et al. 2007). 



M500 (Mq) 

Figure 3. Compajrison of the minimum specific energy required 
to reduce cluster gas mass fractions from the universal value to 
the observed values with estimates of the amount of AGN energy 
available. Solid squares represent the Chandra data of V06 and 
open triangles represent the XMM Newton data of PAP. The three 
dotted lines represent theoretical estimates of i?AGN, defined as 
eMgiic^ where e is the efficiency factor. The thick dashed Une 
shows the amount of energy deposited into the ICM by the AGN 
in Hercules A, which is the most energetic AGN outtjurst known 
(Nulsen et al. 2005). The thin dashed line corresponds to the case 
where a given cluster experiences 10 such Hercules-like outbursts 
over its Ufetime. 



energy by the typical observed mass of gas within r500. In 
particular, we assume Mgas(r5oo) ~ O.I2M500; see Fig. 1.). 
In order to explain the most massive systems, we calculate 
that a minimum energy of ~ 10 keV per particle is required. 
If one adopts an efficiency of e = 0.1, which is approxi- 
mately the efficiency predicted by standard radiatively effi- 
cient accretion disk models (e.g., Shakura & Sunyaev 1973), 
there is potentially just enough energy available in black 
holes distributed throughout rsoo to account for the observed 
gas mass fractions. (We use the term 'potentially' since we 
remind the reader that we have calculated the minimum 
energy required and furthermore have made optimistic as- 
sumptions about the mass of black holes available to heat 
the ICM.) This result agrees quite well with the more de- 
tailed calculations of Bode et al. (2007) (that include, e.g., 
the effects of asphericity and substructure), when one nor- 
malises their estimated required total energy to the same 
total mass of gas Eissumed above. 

However, modelling of AGN-blown X-ray cavities sug- 
gests the typical cluster black hole efficiency is actually much 
lower than 0.1. The most energetic AGN outbursts known, 
in Hercules A (Nulsen et al. 2005) and MS0735.6-f7421 (Mc- 
Namara et al. 2005), have mean powers of « 1.6 — 1.7 x 10*® 
ergs s~^. The typical age of such outbursts is ~ 100 Myr, 
corresponding to a total energy of few times lO''^ ergs or 
a specific energy of a few tenths of a keV per particle (see 
the thick dashed line in Fig. 3). This falls nearly two orders 



3.3 DifFerent cosmological parameters: flm and h 

The observed gas fractions arc proportional to h~^'^, while 
the WMAP constraint is independent of h. Thus, adopting 
a lower value of ft « 0.55 would bring these two results into 
agreement. Indeed, this is why similar analyses by Roussel 
et al. (2000) and Sadat et al. (2005), who adopt h = 0.5, 
find higher gas fractions than we have shown here. However, 
the large body of independent evidence in favour of ft > 0.6 
(e.g. York et al. 2005; Jones et al. 2005; Riess et al. 2005; 
Ngeow & Kanbur 2006) makes this solution seem unlikely. 

We can use the observed baryon content of clusters to 
reverse-engineer the universal total matter density, firn, via: 



6dep(»'500)f26 



fcHSE/gas('"500)[l + Mstar (rsOO ) A^gas (r.soo)] 



(2) 



where &dop('^'50o) is the baryon depletion factor within rsoo 
and fonsE is the observational bias introduced by assuming 
strict HSE (i.e., ignores pressure support due to turbulent 
motions), both of which can be estimated using gasdynamic 
cosmological siirmlations. 

Adopting the WMAP baryon density of Qbh^ = 
0.0223j;Q;QQog (which is in good agreement with the latest 
QSO constraints; see O'Meara et al. 2006), 6dep(^60o) = 
0.95 ± 0.05 (Kravtsov et al. 2005; Ettori et al. 2006; 
Grain et al. 2006), 6hse ~ 0.9 (Nagai et al. 2007), 
Mstar- (rsoo )/Mgas (rsoo) = 0.10 ± 0.05 (see §3.1), and 
/gas(r6oo) = 0.115 ± 0.015 (see Fig. 1 - note that this spans 
both the Chandra and XMM-Newton results at rsoo), we 
find: 

— 0.36^Q 08 

This is larger than the best-fit WMAP 3-year value of Qm = 
0.238 ±0.03, but is in good agreement with several previous 
X-ray analyses (e.g., Allen et al. 2002; Lin et al. 2003; Hicks 
et al. 2006). 

However, our findings differ somewhat from those of Et- 
tori (2003). He found evidence for only a marginal discrep- 
ancy between his cluster baryon fraction estimates (inferred 
from BeppoSAX X-ray data) and the ffrst-year WMAP re- 
sults (Spcrgcl et al. 2003). This is somewhat surprising since 
Ettori (2003) used the data of Ettori et al. (2002), who ac- 
tually find lower gas mass fractions than we have assumed 
above. A likely explanation for this discrepancy is that: (1) 
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he assumed a slightly higher stellar mass fraction; (2) he in- 
voked a large fraction of the cluster baryons being hidden 
at "warm" temperatures of 10^ < T < 10^ K; and (3) he 
extrapolated his results well beyond the maximum observ- 
able radius (~ riooo) to the virial radius assuming a King 
total mass profile for several of his systems. The first two 
are inconsistent with recent, more reliable data (see §3.1). 
Finally, extrapolation of the data out to the virial radius 
is highly uncertain. Since simulations and observations now 
reliably extend out to rsoo, such large extrapolations are no 
longer required. 

In addition to the agreement with several previous X- 
ray studies, our results on fi^ also agree with Sloan Digital 
Sky Survey (SDSS) measurements of the Lyman alpha forest 
power spectrum ( Viel & Haehnelt 2005) , the baryon acoustic 
peak of luminous red galaxies (Eisenstein et al. 2005), and 
the power spectrum of galaxies (Tegmark et al. 2004). We 
also note that combining WMAP 3-year joint constraints 
on ilrnh^ and erg (which are degenerate, see below) with 
weak lensing cosmic shear measurements (Hoekstra et al. 

2006) results in an increased best-fit value for Qm that is in 
good agreement with our results (see Fig. 7 or Spergel et al. 

2007) . Finally, Li et al. (2006) have also recently reported 
that the "low" value of flrn reported by Spergel et al. (2007) 
is in discord with the number of observed strong lensing 
giant arcs, however a fiducial flat model with Qm = 0.3 and 
(78 — 0.9 is able to match the lensing data. 

On the other hand, several other cosmological tests tend 
to support the "low" value of Qm advocated by Spergel et al. 
(2007). These include measurements of the motions within 
the local supercluster (Mohayaee & TuUy 2005), the power 
spectrum of galaxies in the 2dFGRS (Cole et al. 2005), and 
the mass function of galaxy clusters in the HIFLUGCS X- 
ray flux- limited sample (Reiprich 2006). 

Clearly, the relatively large scatter in the reported best 
fit values of Qm warrants further detailed investigation. 
While our own estimate of Qrn = 0.36to'o8 ^ marginally 
inconsistent with the WMAP 3-year constraints (and some 
of the other tests that support the WMAP result), we have 
argued that this discrepancy is a small price to pay for ob- 
taining agreement with the observed cluster gas mass frac- 
tions. Finally, we point out that the discussion carries im- 
plications beyond the precise value of Qm- Since Qm and as 
are strongly degenerate in the CMB data, our results sug- 
gest that as is larger than the best-fit value of 0.74to og de- 
rived from WMAP data alone. It is interesting to note that 
a model with Qm ~ 0.28 and as ~ 0.8 is within ~ 1— sigma 
of both the cluster gas mass fraction data and the WMAP 
joint constraints on these parameters. 



4 CONCLUSIONS 

Recent, good quality observations of massive clusters with 
Chandra and XMM-Newton put strong observational con- 
straints on the gas mass fraction in massive clusters at large 
radius. In many cases the new data allow the fraction to 
be constrained out to rsoo (the radius at which the cluster 
density contrast is 500). Simulations of clusters suggest that 
the cluster gas fraction at this radius should closely reflect 
the average baryon mass fraction in the universe as a whole. 
We find /gas(r5oo) ~ (0.12 ±0.02) (/l/0.72)"^■^ where 



the uncertainty reflects the systematic variations between 
clusters. This is lower than the best fit to the WMAP 3- year 
result of /i, = 0.176 ±0.02. 

We consider whether the discrepancy could be due to 
a large fraction of the cluster gas cooling to form stars and 
cold gas clouds, or whether it could be due to strong non- 
gravitational heating transporting ~ 30% of the cluster gas 
outside rsoo . Observational limits on the stellar and cold gas 
content of clusters appear to rule out the first possibility. 
In order to investigate the second, we compute the energy 
budget required to rearrange the cluster gas. The energy re- 
quired significantly exceeds the plausible energy input from 
black holes, unless their mass accretion history is always as- 
sociated with efficient jet production. 

The most likely explanation is that Qm lies in the range 
0.36lQjg. Such a relatively high Qm lies slightly above the 
68 per cent confidence limits from WMAP, but is consistent 
with the current data at the 2a level. Since the measure- 
ments of Qm and as are highly correlated in the WMAP 
analysis, this means that as is also likely to be higher than 
the formal best-fit value of 0.74. This has important im- 
plications for the abundance of collapsed objects prior to 
re-ionisation (Reed et al. 2007). 
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